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In contrast to the static operations of conventional semiconductor devices, the dynamic confor- 
mational freedom in molecular devices opens up the possibility of using molecules as new types of 
devices such as a molecular conformational switch or for molecular data storage. Bistable molecules, 
with e.g. two stable cis and trans isomeric configurations, could provide, once clamped between two 
electrodes, a switching phenomenon in the nonequilibrium current response. Here, we model molec- 
ular switch junctions formed at silicon contacts and demonstrate the potential of tunable molecular 
switches in electrode/molecule/electrode configurations. Using the non equilibrium Green function 
approach implemented with the density-functional-based tight-binding theory, a series of properties 
such as electron transmissions, I-V characteristics in the different isomer-conformations, and po- 
tential energy surfaces as a function of the reaction coordinates along the trans to cis transition 
were calculated. Furthermore, in order to investigate stability of molecular switches in ambient 
conditions, molecular dynamics (MD) simulations at room temperature were performed and time- 
dependent fluctuations of the conductance along the MD pathways were calculated. Our numerical 
results show that the transmission spectra of the cis isomers are more conductive than trans coun- 
terparts inside the bias window for all two model molecules. The current-voltage characteristics 
consequently show the same trends. Additionally, the calculations of time-dependent transmission 
fluctuations along the MD pathways have shown that the transmission in cis isomers is always sig- 
nificantly larger than that of trans counterparts showing that molecular switches can be expected 
to work as robust molecular switching components. 



I. INTRODUCTION 

Since the successful measurement of current flow 
through individual molecules in the past few decades, 
various kinds of molecular devices have been demon- 
strated by many groups, showing interesting proper- 
ties such as negative differential resistance 1 and mem- 
ory effects. 2 The realization of electronic function at the 
molecular scale, visionally proposed by Aviram and Rat- 
ner already in 1974^ is gradually finding experimen- 
tal support due to the development of scanning probe 
technologies^ Nevertheless, many challenges are still 
ahead. The conductance of a molecule connected be- 
tween contacts changes dramatically depending on many 
parameters such as the coupling strength between the 
molecule and of the contacts, the electronic structure of 
the molecule and contacts, and molecular conformational 
changes. Therefore, it is highly desirable to understand 
the relationship between those influential parameters and 
transport properties for the design of controllable molec- 
ular devices. 

In contrast to the static operations of conventional 
semiconductor devices, the dynamic conformational free- 
dom in molecular devices opens the possibility of using in- 
dividual molecules as a new type of device such as a con- 
formational molecular switch or for molecular data stor- 
age. The ability to control the conductance at the single 
molecule scale will have a great impact in nanoscale elec- 
tronics as well as other fields such as chemical sensors^ 
or medical diagnostics. 6 Photochromic organic molecules 
having bistable characteristics are attractive materials 
for the design of the molecular switches since a great 
variety of properties such as their molecular conforma- 



tions and their electronic structures can be reversibly al- 
tered from one state to another by external stimulations 
such as light irradiation and current-pulses. Azoben- 
zene, diarylethene and their derivatives are representa- 
tive photochromic molecules. Recently there have been 
many theoretical and experimental studies on molecu- 
lar switches i Theoretical studies have explored the use 
of azobenzene as a molecular switch when connected by 
Au electrodes^ or carbon nanotubes (CNTs)^ Promising 
results were obtained, as a significant change in conduc- 
tance was seen between the two conformations. Addi- 
tionally, the study of diarylethene derivatives connected 
to Au contacts also have presented large change in con- 
ductance between open or closed formsJ^ 

Several experimental groups have addressed the fabri- 
cation of molecular switches consisting of molecules with 
bistable characteristics sandwiched between contacts. 
Dulic et at. have created diarylethene derivative-based 
molecular switch junctions connected between gold con- 
tacts by means of the mechanically controlled break junc- 
tion (MCBJ) technique and measured the current flow 
through the junction before and after UV irradiation #AI 
The current flow through the molecule was suddenly in- 
creased by the UV irradiation. Additionally, the current 
flow did not decay after stopping the irradiation. Al- 
though the reversible switching has not been observed, 
conformational switching and the control of current flow 
through the single molecular device have been realized. 

Lortscher et at. have fabricated molecular junctions 
where single bipyridyl-dinitro oligophenylene dithiole 
molecules have been contacted between gold contacts 
with the MCBJ technique^ They have demonstrated re- 
versible molecular switching by applying current pulses 
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across the gold contacts. Additionally, Henzl et ali^ and 
Choi et ali^ have reported reversible cis-trans isomer- 
ization of a single azobenzene derivative at the Au(lll) 
surfaces induced by the inelastic excitation of vibrational 
states by an STM tip. The experiment by Alemani et al. 
have reported mechanical reversible isomerization of a 
single azobenzene derivative on the Au(lll) induced by 
an electric field applied from an STM tip not in contact 
with the molecule. 15 

In this article, we explore the potential of azobenzene- 
based molecules as reversibly tunable molecular switches 
in electrode/molecule/electrode configurations. We have 
modeled two kinds of molecular switch junctions formed 
at silicon contacts. In order to improve the molecule- 
contact coupling the molecules are covalently attached 
to Si contacts. We present a detailed analysis of the en- 
ergetics, potential energy surface (PES) as a function of 
the reaction coordinate along the trans to cis transition, 
the stability along MD pathways, and transport prop- 
erties in the two conformations using NEGF approach 
based on DFTB method. 

This paper is organized as follows. In Sec. II, we 
present computational details which are used in our 
study. In Sec. Ill, we show the molecular structures, 
electronic structures, and PESs of molecular switches in 
gas phase and connected between Si contacts. In Sec. IV, 
we show the transport properties and the stability of the 
molecules between Si contacts. We conclude in Sec. V. 



DFTB method to include the NEGF formalism com- 
bined with the Landauer formalism. The gDFTB code 
is based on a tight-binding representation of the wave 
function and a density-functional formulation of the 
Coulomb and exchange correlation potentials. The 
conductance at low or zero bias (equilibrium) is ob- 
tained from electronic transmission probability T as 
G = G T(E,V = 0), where G = 2e 2 /h is the quan- 
tum conductance, e is the electron charge, and h is 
the Plank constant. The electronic transmission prob- 
ability is related to the retarded/advanced Green func- 
tion of molecular region G R ' A {E) and self-energies of 
the electrodes E L / R (£?) by Fisher and Lee relation^ 3 - 
T(E,V) = 4Ti[G R (E)lm^ L (E)G A (E)lmJ: R (E)]. We 
emphasize that when calculating the current I(V) = 
(2e/h)JdE[f L (E) - f R (E)]T(E,V), the full voltage- 
dependent T(E, V) is calculated through self-consistent 
Poisson-NEGF cycle in gDFTB. Here, f L/R (E) are the 
Fermi functions of the left/right electrode. The minimal 
basis set used in gDFTB could be responsible for inaccu- 
racies in the calculation of the current. This may result 
in problems describing transport mediated via the tails of 
the wave functions. This problem could be solved by us- 
ing an ab initio code with a richer basis set, which would 
be dramatically increase the computational requirements 
without adding any qualitative new results to our treat- 
ment. 



II. METHOD 



The relaxation of the molecular structures, total en- 
ergy calculations, and electronic structure calculations 
are based on a second-order expansion of the Kohn-Sham 
total energy in density-functional theory with respect to 
charge density fluctuations^ The zero-th order approach 
is equivalent to a common standard non-self-consistent 
tight-binding scheme, while at second order a transpar- 
ent, parameter- free, and readily calculable expression for 
generalized Hamiltonian matrix elements can be derived. 
The resulting DFTB method provides accurate results 
comparable to the results obtained with ab-initio DFT 
methods or first principle calculations while requiring less 
computational resources. These are modified by self- 
consistent redistribution of Mulliken charges. The com- 
putations in DFTB+ use a minimal basis set, two cen- 
ter approximation, neglect of three centers and crystal 
field. The atomic pair potentials are computed at a pre- 
computation stage using a DFT code. The Slater-Koster 
integrals are stored (as a function of the inter-atomic dis- 
tance) in look-up tables which are used during the actual 
simulation. DFTB+ have been extensively validated in 
the past years by many applications related to e.g. semi- 
conductor surface reconstructions j 17 i 18 i 19 reaction of or- 
ganic molecules ^ 18 ! 20 and biomolecules i 18 i 21 

Conductance calculations were done via gDFTB 
code* 2 ^ The gDFTB code is the extension of the 



III. PROPERTIES OF MOLECULES IN GAS 
PHASE AND BETWEEN SI CONTACTS 



Figure 1(a) shows the two model systems considered 
in this work. Two different kinds of pairs of linkers are 
attached to the azobenzene molecules. We refer to these 
two model systems named as system A and system B. 
In each system, two linkers covalently connecting the 
molecule to silicon contacts are substituted at meta and 
ortho positions of the phenyl rings. The large movement 
of one phenyl ring involved with the cis-trans isomeriza- 
tion is unfavorable for the stable operation of molecular 
switching since such large movement prevents stable cou- 
pling between the molecule and Si contacts. However, 
via the linkers at the meta and ortho positions unfavor- 
able displacement at molecule-contact interface can be 
reduced as shown in Fig. 1(b). Additionally, the change 
of molecular lengths from one isomeric state to the other 
is small (see Table I) so that the molecules do not tilt 
through the cis-trans isomerization. The atoms at con- 
tacts do not move throughout, during the cis-trans iso- 
merization. Thus, the choice of these positions of substi- 
tutions is preferable for the stable operation of molecular 
switching. 
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System A (0-0) System B (N-N) 

cis 10.208 A 11.383 A 

trans 11.000 A 10.685 A 

TABLE I: Molecular length of two molecular switches in gas 
phase: The molecular lengths are defined by the atomic pairs 
in parentheses which make bonds with silicon surfaces at top 
and bottom. These values are obtained after structure op- 
timization by DFTB+ calculations. For the detail see the 
text. 



(a) (b) 




System A System B 



FIG. 1: (a) Two azobenzene derivatives modeled in this work 
and (b) schematic of azobenzene-based molecular switch con- 
nected to Si contacts via linkers. The two linkers are sub- 
stituted at meta and ortho positions of phenyl rings. These 
positions of substitutions reduce the movement of atoms at 
contact involved with cis-trans isomerization. 



A. Isolated molecules 

For the determination of stable configuration of the 
molecules between Si contacts, we first performed a ge- 
ometry optimization for the isolated molecules by means 
of the conjugate-gradient technique. The geometry opti- 
mization was carried out until the absolute value of the 
interatomic force reduced to less than 10 -4 atomic units. 
We used the Slater-Koster parameters developed by El- 
stner et at. for the C, H, N, O atoms^ and by Szuecs et 
al. for the Si atom as basis setsJ^ Figure 2(a) and 2(c) 
show the optimized structures, their frontier orbit als as 
well as the molecular energy levels of system A and B, 
respectively. The trans isomers are characterized by a 
planar structure while in each cis isomer two phenyl rings 
are twisted each other of 12.9 degree in system A and B. 

Figure 3(a) shows the PES of total energy as a function 
of the reaction coordinates along the trans to cis tran- 
sition for the two model systems in the gas phase. The 
total energy calculations have shown that trans isomers 
in system A and B are energetically more stable than cis 
counterparts by 0.193 eV and 0.277 eV, respectively. The 
trans to cis potential energy barriers for system A and B 
are 1.060 eV and 1.433 eV, respectively. These features 



are similar to the azobenzene molecule in the gas phase 
where its trans isomer is experimentally reported to be 
more stable than the cis counterpart by 0.6 eV and the 
potential energy barrier from the trans to the cis coun- 
terpart is 1.0 eV*2i 



B. Molecules between Si contacts 

Let us consider the model systems connected between 
Si contacts. Hereafter we refer to these coupled sys- 
tems as "extended molecules". The complexity of the 
systems leads to a very complicated PES with many 
metastable states; for simplicity we only consider the 
vertical arrangements of the molecules with respect to 
the hydrogen-passivated Si(lll) surface. The structure- 
optimized molecules in the gas phase are covalently con- 
nected to hydrogen-passivated Si(lll) surfaces and re- 
laxed again using DFTB+ with periodic boundary con- 
ditions. The spacing between lateral organic molecules is 
set large enough to exclude lateral interactions between 
azobenzene molecules. In each unit cell, the top/bottom 
silicon contacts are comprised of 54 Si atoms forming 
three silicon layers and 8 hydrogen atoms covering the 
Si(lll) surface. In order to include the reconstruction of 
the silicon atoms at the surfaces two silicon layers from 
the surface are relaxed while the third layers of silicon 
contacts are fixed, which is required for the subsequent 
conductance calculations. The Si-Si bond lengths in the 
fixed layers are set as 2.35 A, which is experimental value 
for bulk Si crystal. Figure 4 shows the total energy of the 
extended molecules as a function of the separation of two 
Si surfaces. Figure 5 shows the optimized unit cells of two 
model systems at the separation of Si surfaces where the 
total energies of trans isomers between Si contacts are 
their local minimums which are shown in vertical dot- 
ted lines in Fig. 4. The PES calculations of extended 
molecules and subsequent conductance calculations are 
carried out at these separations. 





System A 


System B 


cis 


-4.259 eV 


-4.255 eV 


trans 


-4.259 eV 


-4.253 eV 



TABLE II: Fermi energies of two molecular switching junc- 
tions. 

The PESs of total energies in the extended molecules 
are shown in Fig. 3(b). The potential energy barriers of 
extended molecules along the reaction coordinate from 
trans to cis isomers in system A and B are 2.001 eV 
and 1.675 eV, respectively. The potential energy barriers 
are increased compared with those in the gas phase (See 
Fig. 3(a)). This is due to the reduction of the conforma- 
tional degree of freedom of molecules confined between 
Si contacts. 



System A 




-5.30 eV -5.31 eV -5.23 eV _ 5 , 32 e v 



System B 




HOMO-7 HOMO-8 HOMO-7 HOMO-8 

-5. 1 5 eV -5.25 eV -5.1 1 eV -5.22 eV 



FIG. 2: Energy level diagram and frontier orbitals of system A in (a) gas phase and (b) between Si contacts. Energy level 
diagram and frontier orbitals of system B in (c) gas phase and (d) between Si contacts. 
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FIG. 3: Potential energy surface of both systems A and B (a) 
in gas phase and (b) between Si contacts. The total energies 
of both systems are plotted as a function of the dihedral angle 
between the two bonds connecting nitrogen atom and phenyl 
rings. The curves interpolates points, which are indicated by 
the crosses. 



FIG. 4: Total energies as a function of distance between two 
silicon surfaces for two molecular switch junctions. The bro- 
ken lines indicate the distances where the molecules in trans 
conformation are energetically most stable. All transmission 
calculations in this article are carried out with these distances. 



5 



(a) System A 



(b) System B 




FIG. 5: Optimized structure (unit cell) of two molecular 
switch junctions modeled in this work. cis (trans) con- 
formations are shown in red (green). For comparison, the 
two isomer conformations are superimposed. Optimizations 
are carried out using the conjugated-gradient technique with 
DFTB+ under periodic boundary conditions. The silicon lay- 
ers at the bottom and top are fixed. 
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IV. TRANSPORT PROPERTIES 

Once the relaxed geometries have been determined, we 
first focus on the electronic transmission probabilities of 
molecular switch junctions at equilibrium for cis and 
trans isomer-conformations. The structure-optimized 
unit cell of the extended molecules is put between semi- 
infinite Si bulk layers and then transmission calculations 
are carried out. Figure 6 shows the electronic transmis- 
sion probabilities at equilibrium. The Fermi energies are 
near the edge of valence band (for the Fermi energies, 
see Table II). Thus, the conductance in these model sys- 
tems is driven by hole transport under small bias. Our 
results show that the cis isomers are more conductive 
than the trans counterparts around the Fermi energies 
in both of the model systems. The low conductance be- 
tween the resonant peaks indicates that the transport 
can only take place through tunneling. The conductance 
around the Fermi energy in cis isomers is higher than 
that of trans counterparts by a few orders of magnitude. 
This suggests that switching behavior may be expected 
at low bias. The low transmissions around the Fermi en- 
ergies are due to the weak coupling between the molecules 
and the Si contacts giving rise to narrow resonant peaks. 
As shown in the DOS profiles of both system A and B in 
Fig. 7, the narrow resonant peaks due to the HOMO and 
LUMO of isolated molecules are located at about —5.0 
eV and —3.0 eV, respectively. The Fermi energies lie in 
the center of gaps between these narrow resonant peaks 
due to HOMO and LUMO of isolated molecules. Hence, 
the DOS on the molecules around the Fermi energies are 
very low. 



FIG. 6: Electronic transmission probabilities for system A 
and B at equilibrium. Broken lines show the Fermi energies. 
In (a), Fermi energies for two isomers are same. Thus, only 
one broken line is shown. For the detail of Fermi energies, see 
the Table II. 



Previous theoretical studies of azobenzene derivatives 
connected between Au 8 or CNT— electrodes have re- 
ported that cis isomers are less conductive than their 
trans counterparts since the ^-conjugation in the trans 
isomers along the molecular framework will be broken by 
the torsion involved with cis-trans isomerizations. How- 
ever, the same discussion does not hold true in our case 
since the HOMO of the trans isomers does not consist 
of p z orbitals perpendicular to the molecular framework, 
but non-binding orbitals mainly localized around nitro- 
gen atoms^ as shown in Fig. 2(a) and Fig. 2(c). The 
eigenstates consisting of p z orbitals of the trans isomers 
lie in the HOMO-3 levels. Thus, these states do not 
contribute to the transport at low bias. This counter- 
intuitive result (I (cis) > I (trans)) can also be under- 
stood by the analysis of resulting mixture states of the 
extended molecular systems as follows. 

Figure 6 shows that the transmission curves for cis 
isomers near the Fermi energies are higher than those of 
trans counterparts, implying higher current through cis 
isomers than the trans counterparts at low bias. The 
difference between the transmission curves of the cis and 
trans isomers near the Fermi energies comes from the dif- 
ference between the height of the resonant peaks of the 
cis and trans isomers around —5.0 eV. The height of the 
transmission resonant peaks in the two different isomer- 
conformations can be interpreted from the molecular or- 
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FIG. 7: TDOS for system A and B at equilibrium. TDOS are shown in red curves. The contribution from Si surfaces to 
TDOS is shown in green curves. Broken lines show the Fermi energies. 



bitals of the extended molecules. The HOMOs of the 
cis isomers interact with surface states of the Si contacts 
and their orbit als hybridize. The resulting mixed states, 
corresponding to the HOMO-8 of systems A in Fig. 2(b) 
and of system B in Fig. 2(b), have a non- negligible orbital 
distribution bridging across the Si states. These states 
can contribute to transport as shown in the resonance 
peaks (in red) of system A at —5.31 eV in Fig. 6(a), of 
system B at —5.17 eV in Fig. 6(b). 

On the other hand, the HOMOs of the trans isomers 
have neither orbital amplitude around the linkers nor 
interact strongly with the surface states of the Si con- 
tacts. This leads to less orbital hybridization as can be 
seen for HOMO-8 of Figures 2(b) and (d). These states 
do not contribute to transport as shown by the sharper 
and smaller resonant peaks (in green) for system A at 
—5.22 eV in Fig. 6(a) and for system B at —5.10 eV 
in Fig. 6(b). Despite the delocalized molecular orbitals 
along the molecular framework, the LUMOs of both the 
trans and cis isomers in all model systems do not interact 
with the surface state of the Si surface. This is because 
the energetic difference between LUMO of both cis and 
trans isomers and conduction state of the Si surface (the 
conduction band is above —2.70 eV and valence band is 
below —4.51 eV) is large. Additionally, there is no state 
available in the silicon contacts on the energy around 



LUMOs of both the trans and cis isomers in both of the 
model systems. Thus, these states lead to much lower 
peaks around —3.0 eV in transmission spectra as shown 
in Fig. 6. 

We have to pay attention not only to the 7r-conjugation 
but also to the interaction between the molecule and the 
contacts. These results support the conclusion that the 
electronic structure of both molecule and contact are cru- 
cial in determining the transport properties at low bias. 

In order to confirm the possible switching behavior al- 
ready suggested by the analysis of transmission proba- 
bilities at equilibrium, we have calculated the current- 
voltage (I-V) characteristics for the two different isomer- 
conformations of the two model systems. The I-V char- 
acteristics for cis and trans conformations for the two 
model systems are shown in Fig. 8. The charge trans- 
port in trans isomers is strongly suppressed compared to 
that in the cis counterparts for both of the two model 
systems. If the bias voltage is less than 1.0 V, the reso- 
nant peaks do not enter Fermi bias window so that the 
I-V curves are smooth within the bias window. There- 
fore, we can expect constant on-off ratios as a response 
to the change or fluctuation of bias voltage, which are 
preferable for the robust control of the molecular switch- 
ing. The bottom panels in Fig. 8(a) and (b) show the 
on-off ratios (the ratio of current in cis isomer to trans 
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FIG. 8: I-V characteristics of molecular switch junctions in 
two different conformations: The bottom panels show the on- 
off ratios (the ratio of current between cis and trans isomer). 
In all cases, the cis conformations (red) are conductive than 
trans counterparts (green). 



counterpart). The on-off ratio reaches up to 647 for sys- 
tem B at 0.2 V (bottom panel of Fig. 8(b)). Although the 
on-off ratios decrease with higher bias voltage, all molec- 
ular switches recorded on-off ratios over 10 within 1.0 V. 
Higher on-off ratios for lower bias voltage are advanta- 
geous both for power saving and for the suppression of 
unexpected breakdown triggered at higher bias voltage. 
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FIG. 9: Time dependent fluctuations of transmission spectra 
along MD pathways at equilibrium for two molecular switch 
junctions in Fig. 1. In each graph, 100 transmission spectra 
are superimposed. The time interval of 100 curves is 100 fs. 



A. Time dependent fluctuation of conductance 

For the reliable operation of molecular switching one 
needs to investigate the stability of the model systems 
between Si contacts and the fluctuation of the current 
through them in ambient conditions. For this purpose, 
we have performed MD simulations of the model sys- 



tems between Si contacts and calculated the conduc- 
tance of the two different isomer-conformations along the 
MD pathways. The calculation setting was as follows. 
We used the Verlet algorithm for the MD simulation in 
DFTB+ with periodic boundary conditions. We used the 
same SK parameters used for the structure optimization. 
The MD duration is set for 10 ps, i.e., 10,000 steps at in- 
terval of 1.0 fs. The nuclear and electronic temperatures 
are set as 300 K. After performing MD simulations the 
coordinates in every 100 steps are extracted from MD 
trajectories. Finally we have repeated conductance cal- 
culations for collected 100 coordinates using the gDFTB 
code both with and without applied bias voltage. We 
then analyzed the fluctuation of the transmission proba- 
bilities and I-V curves along the MD pathways. 
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FIG. 10: The dispersion of currents along MD pathways as a 
function of bias voltage. In each bias voltage 100 values of the 
current in cis isomers (trans isomers) along the MD pathways 
are plotted in red (green) crosses. 

The transmission probabilities along the MD pathways 
are plotted in Fig. 9. In each isomer 100 transmission 
spectra are plotted. We can see that the main features 
of the transmission spectra in Fig. 6 remain unchanged. 
Hence, we can expect that the fluctuation of current flow 
under ambient conditions is small. The dispersion of the 
current along the MD pathways as a function of bias volt- 
age is plotted in Fig. 10. The dispersion of the current in 
the two isomers is distinguishable, with the cis isomers 
always being more conductive than the trans counter- 
parts. The currents in the two isomers as a function 
of time under bias (0.2 V) in Fig. 11 also support this 
conclusion. The on-off ratio is over 15 (the minimum 
value is 15.5 at 4.7 ps in Fig. 11(b)) for all model sys- 
tems. Therefore, the molecular switch junctions modeled 
on this study may be expected to work as switching com- 
ponents in ambient conditions. 

Finally we would like to comment on the possi- 
ble experimental realization of molecular conformational 
switching device consisting of azobenzene derivative and 
Si contacts. The attachment of aromatic compounds to 
the Si surface will be possible using the methodology de- 
veloped by Murata et al 26 ^ 27 The realization of two ter- 
minal Si/molecule/Si devices by covering of the function- 
alized Si surfaces from the top with a Si contact might 
present a formidable challenge due to the intrinsic asym- 
metry in molecular junctions. However, first attempts 
to such a difficult task proved already successfully^ 
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FIG. 11: Time dependent fluctuation of current (top) and 
on-off ratio (bottom) along MD pathways in molecular switch 
junctions. 



For further understanding and control of the molecu- 
lar switching, systematic studies investigating the effect 
of contacts, 29 internal scattering effect ^tt surface rough- 
ness, and impurities^ will be required. The spacing be- 
tween two Si surfaces will also be an important factor 
since it will influence the PES, molecular and electronic 
structure, and resulting conductivity of the molecule em- 
bedded between Si contacts. Concerning the design of 
azobenzene derivatives, further investigation is required 
to determine which linkers and functional groups lower 
the PES barrier for cis-trans isomerization. One pos- 
sibility to achieve it is to increase the flexibility of the 
molecule by using long linkers which should not, however, 
increase the fluctuation of current due to the flexibility 
of the molecules. 



V. CONCLUSION 

In summary, we have modeled two molecular switch 
junctions consisting of azobenzene derivatives and sil- 
icon contacts and explored their potential as molec- 
ular switches by the analysis of the energetics, PES, 



and transport properties in the two different isomer- 
confor mat ions: cis and trans. Modeling of the molecu- 
lar switches demonstrated that the attachment of linkers 
at the meta and ortho positions could reduce the steric 
repulsion between the molecules and silicon contacts dur- 
ing cis-trans isomerization to a value which is reasonable 
for the robust operation of molecular switching. Our nu- 
merical results show that cis isomers are more conduc- 
tive than the trans counterparts inside the bias window 
for both model systems. I-V characteristics also show 
the same trends. Although the 7r-conjugation through 
the molecular framework is reduced due to the torsion 
of two phenyl rings, the cis isomers interact much more 
strongly with the surface states of the silicon contacts 
than their trans counterparts, resulting in larger conduc- 
tivities. Additionally, the calculation of currents through 
the two model systems along the MD pathways in ambi- 
ent conditions revealed that the currents in cis isomers 
are always significantly larger than trans counterparts 
leading to high on-off ratio. Importantly, the disper- 
sion of the current as a function of time is small enough 
to distinguish two different isomer-conformations. Thus, 
azobenzene derivatives between Si contacts are promis- 
ing for applications as robust molecular switching com- 
ponents or single molecular data storage devices. 
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